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Abstract 

Factorial and cumulant moments in a simple cascade models are considered. 
Their characteristic features are shown to be similar to those observed in QCD 

jets. 



1 Introduction 



The characteristics of multiplicity distributions in multiparticle dynamics has re- 
cently gained much attention due to a QCD analysis of generating functional 
of these distributions and their moments. A very clear review of this subject 
can be found in Q, where one can also find the relevant references, see also a 
recent review 0. The main focus in the above studies was initially a striking 
QCD prediction on the oscillatory behavior of cumulants as functions of their 
rank. Actually it turned out to be more convenient to analyze a ratio of cumu- 
lants to factorial moments. The experimental data from many reactions (see @) 
confirmed this prediction. The situation seemed to be all the more interesting, 
that conventional multiplicity distributions such as negative binomial one, did 
not predict such oscillations. However, a subsequent analysis || has shown, that 
other more phenomenological models such as dual parton model [|J or quark- 
gluon string model || are also explaining the data. This brings us to a natural 
question. What is the real (basic) origin of the oscillations of the cumulants? 
To answer this question we construct a simple probabilistic cascade and look at 
the properties of the cumulants. We find exactly the same oscillation pattern as 
predicted theoretically and observed experimentally, thus confirming one of the 
conclusions of that the basic origin of these oscillations is a cascading origin 
of the underlying process of particle formation. 



Let us consider a model of cascade with discrete steps in time. Initially (at time 
t = 0) we have only one "particle". When time changes from t to t + 1 any 
particle existing at a time t splits with probability p (a parameter of the model) 
into 2 particles or does not split with probability 1 — p. 

Now P n (t) is a probability that on t step n particles were produced. The "first 
split" equation reads 



P n (t) = (l-p)P n (t-i)+p y E P ni(t-l)Pn- ni (t-l) , Pn(0) = 5 nl (1) 



Let us now consider a generating function related to the probabilities P n (t): 



2 The model. 



n 



ni=0 



oo 



n 



(2) 



n=0 



From ([I]) we have 

G(t, z) = (1 - p)G(t - 1, z) + P G 2 (t - 1, z) 



(3) 



The boundary conditions for this relation is 

1. G(t, 0) = 1 - the total probability conservation; 

2. G(0, z) = 1 + z - corresponding to P„(0) = S n i 
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3 Mean multiplicity and moments. 

By definition a mean multiplicity n(t) is related to a generating function by 

dG{t, z 



n{t) = Y j nP n (t) 

n=0 OZ 



(4) 

2=0 



It gives for (|3|) 

n (t) = (1 - p)n{t - 1) + 2pn(t - 1) = (1 + - 1) = (1 + pf. (5) 

Last equality follows from G(t, 0) = 1 (see above) and n(0) = 1. 

Factorial moments are introduced as coefficients of F q (t) in the expansion 

oo q 

G(t,z) = l + Y,-FMn q (t). (6) 

9=1 I' 



After substitution into (|^) we find 

I 

k=0 



1 + pfF q {t) = (1 - p)F q (t - 1) + p ]T CjF A (t - - 1) , g > 1 (7) 



It is easy to demonstrate that the limit ^lim F q (t) = F q exists. Therefore from 
© we get 



p 1 p 



5-1 



^ — y C k n F k F Q _ k , q > 1 Fi = 1 (8) 

Now it is possible to treat asymptotic factorial moments as factorial moments 
of " asymptotic" generating function. 

oo p 

G ( a ) = 1 + E "f a9 - ( 9 ) 

The equation on G(a) now reads 

G((l+p)a) = (l-p)G(a)+ P G 2 (a) (10) 
with the initial conditions G(0) = 1 and G'(0) = 1. 



4 Cumulants and 

Cumulants are introduced as coefficients K q (t) in the expansion 

oo q 

\nG{t,z) = Y J -7K q (t)n\t). (11) 

q=l I' 
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The cumulant and factorial moments are actually not independent. The relation 
between the two can be traced from the formula G z (\nG)' z = G' z . From @ and 
(|TTD we obtain 

K q = F q - J2 Cl_ x K h F^ k (12) 
fc=i 

It turns out convenient to consider the ratios H q = K q /F q obtained from (|8|) 
and (0). The resulting H q for some values of p are shown in Fig. 1. We see, 
that the ratio of cumulant to factorial moments H q in our simple cascade model 
provides the same oscillation pattern that follows from a number of theoretical 
models and is observed experimentally. How good is a quantitative predictions of 
the model when compared to the experimental data is questionable. We restrict 
ourselves to illustrate this situation at Fig. 2. 

It is instructive to look more attentively at two extreme limits of our cascade 
model. 

1. p -»• 

Keeping the leading power of p in (llOD we get 

paG\a) = pG(a)(G(a) - 1) (13) 
The solution of this equation that satisfies boundary conditions is 

G(a) = - J- (14) 
1 — a 

so F^ = q\, = (q — 1)! and = 1/q. It is worth noting that this result 
is valid only for pq -C 1 (not for p<l). 

2. p = 1 

Although in the context of our study this limit looks quite peculiar, let us proceed 
and perform an explicit calculation, which gives 

G{2a) = G 2 {a). (15) 

The obvious solution of this equation is 

G(a) = e a (16) 

so F^ = 1, KfV = 5 q \ and = 5 q \. Such moments are produced by Poisson 
distribution. 



5 The continuous limit. 

In the context of our model investigations it is useful to determine continuous 
limit of (Q). This allows to calculate not only asymptotic shape of generating 
function and moments but also their dependence on mean multiplicity n. 
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It is convenient to give (|3|) the form 

G(t, z) - G(t -l,z)= pG[t - 1, z){G{t -l,z)- 1). (17) 

If p is considered to be a small one can replace the difference in l.h.s. of this 
equation by derivative (the change at one step becomes infinitesimal) 

^^l=p G (t,z)(G(t,z)-l) (18) 

(we have also changed t — ltot). In this limit the equation on mean multiplicity 
is (see |]) 

d ^ = Mt) (19) 

that is 

n(t) = e pt n(0) = e pt (20) 
according to the initial condition P„(0) = 8 n \. Solution of (|T8|) is clear 

G(M) = (i-=e>£)", (2i) 

where initial conditions and (^) have been used. From this we obtain factorial 
moments and cumulants 

i^C*) = ffl f 1 " -LV (22) 



K q (t) = (q-1)\ 1--7TS -37777 ( 2 3) 



nit), 








n{t)j 


' 71' (t) 


1 





SO 

h ^ = K 1 - 4j i 1 - rt^j • (24) 

In the limit n — > oo this result becomes the same as in (0). 

Here it is appropriate to discuss some relations between our model (in contin- 
uous limit) and a PB (pure birth) one considered in ||. That model is formulated 
as an equation on probability to obtain n particles at the moment of (continuous!) 
time t 

dP(n t) 

— 7T— - — - = X(n - l)P(n - 1, t) - XnP(n, t). (25) 

For this equation on P's the equation on generating function (which we denote 
by G) is 

dG x ,dG 

■W = Xz(1 + z) J7 (26) 
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Mean multiplicity n in this case satisfies 

\n(t) = e M h(0). (27) 




dt 

It is easy to see that the equation ( Effi) has a solution 

-/(?^). (-) 

Here / is a function determined by initial (t = 0) distribution. Universal con- 
ditions on / are only /(0) = 1 (from G(t,z = 0) = 1) and /'(0) = 1 (from 
dG(t, z = 0)/dz = n(t)). So in this model the asymptotic (n — > oo) behavior of 
factorial moments and cumulants is determined by initial distribution while the 
dependence on n is described by differential equations (P5| ) and (|26|). 

It is clear that in the continuous limit our model turns into PB if one specifies 
some initial conditions to (13). 



6 Conclusions 

From the study of the simple cascade model of particle multiplication we see, that 
it correctly reproduces the oscillation pattern predicted by a number of theoretical 
schemes and observed experimentally. We conclude that the basic origin of such 
oscillations is a cascade structure of underlying dynamic picture. This cascade 
could in principle originate both from the quark-gluon phase of the evolution of 
hadronic system and from the nonperturbative hadronization stage. 

We are grateful to I.M. Dremin and I.V. Andreev for useful discussions and 
suggestions. This study was supported by Russian Foundation for Basic Research 
under Grant 96-02-16210a. 



References 

[1] I.M. Dremin, Usp. Fiz. Nauk 164 (1994), 785 



[2] I.M. Dremin, Multiplicity Moments in QCD and Experiment, |hcp- 
| ph/9604245] 1996 

[3] A. Capella, I.M. Dremin, V.A. Nechitailo and J. Tran Thanh Van, Moment 
Analysis of Multiplicity Distributions, |iep-ph/9604247| , 1996, Zeit. Phys. C, 
to be published 

[4] A. Capella, J. Tran Thanh Van, Z. Phys. C23 (1984), 165; A. Capella et al., 
Phys. Rev. D35 (1987), 2921 



6 



[5] A.B. Kaidalov, Sov. J. Nucl. Phys. 45 (1987), 902 

[6] M. Biyajima et al. Phys. Rev. Lett. B237 (1990), 563 



7 




Fig. 1 Cumulants to factorial moments ratio in the model for 3 values of p. 




Fig. 2 



